PRESCRIBING CURVATURES ON THREE DIMENSIONAL 
RIEMANNIAN MANIFOLDS WITH BOUNDARIES 



LEI ZHANG 

Abstract. Let (Af, g) be a complete three dimensional Riemannian 
manifold with boundary dM. Given smooth functions K{x) > and 
c{x) defined on AI and dM, respectively, it is natural to ask whether 
there exist metrics conformal to g so that under these new metrics, K 
is the scalar curvature and c is the boundary mean curvature. All such 
metrics can be described by a prescribing curvature equation with a 
boundary condition. With suitable assumptions on K,c and (M, g) we 
show that all the solutions of the equation can only blow up at finite 
points over each compact subset of M, some of them may appear on dM. 
We describe the asymptotic behavior of the blowup solutions around 
each blowup point and derive an energy estimate as a consequence. 



1. Introduction 

In geometric analysis two well known problems are closely related, the 
Nirenberg problem (or the Kazdan- Warner problem) and the Yamabe prob- 
lem. The Nirenberg problem asks what function K{x) on is the scalar 
curvature of a metric g on S" conformal to the standard metric gQ. From an- 
other point of view, on any given compact Riemannian manifold ( Af, g) with- 
out boundary, the Yamabe problem, which was solved through the works of 
Yamabe [36], Trudinger[35], Aubin[3] and Schoen [32], concerns whether it 
is possible to deform g conformally to get a new metric with constant scalar 
curvature. Similar questions can also be asked on general Riemannian mani- 
folds with boundaries: Let (M, g) be a Riemannian manifold with boundary, 
let K G C^(M) and c G C^{dM) be functions defined on M and dM, 
respectively, then is it possible to deform g conformally to another metric 
gi so that K is the scalar curvature and c is the boundary mean curvature 

4 

under gi? By writing gi as gi = u^-'^g, K and c are related to the scalar 
curvature R„ and the mean curvature /i„ under metric g by 
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where Vg is the unit normal vector pointing to the outside of dM. is the 
Laplace-Beltrami operator, which can be written as A^ = -^di{^g^^ dj) in 
local coordinates. 

If K and c are constants and {M,g) is compact, the existence of a met- 
ric with constant scalar curvature and boundary mean curvature is always 
referred to as the boundary Yamabe problem (BYP). Many cases of the 
BYP have been solved by Escobar [III [THl HI] , Han-Li[22] and Marques [28]. 
But unlike the completely solved Yamabe problem, some cases of the BYP 
are still open. In general, the boundary terms in (jl.ip make the nature 
of equation (jl.lh very different from its counterpart without the boundary 
condition. Some difficulties created by the boundary terms are still not com- 
pletely understood. In this article, we focus on three dimensional Riemann- 
ian manifolds with boundaries and consider the corresponding prescribing 
curvature equations defined on these manifolds. 

Let (M, g) be a smooth three dimensional complete Riemannian manifold 
with boundary dM. Suppose K{x) > and c(x) are functions defined on 
M and dM, respectively. If 171 = u'^g (u > smooth) is a metric conformal 
to g that takes K as the scalar curvature and c as the mean curvature on 
dM, one can write the equation as 



The purpose of this article is to understand the bubbling phenomena of (jl.2p 
under natural assumptions on K and c. Without the boundary condition, 
the bubbling phenomenon of equation (jl.2p and its high dimensional variants 
have been studied extensively under various assumptions. The reader is 
referred to the following incomplete list and the references therein [H O [6l 



EllIIlliaiiaiiaillllSOlllIllSaiMlIlHlIMll^ However, 



much less work can be found to address the case with boundary conditions. 
In this article we describe the blowup phenomenon for (jl.2p under weak 
assumptions on K and c. 

To state our main result, we first remark that M may not be compact, so 
our result concerns the blowup phenomenon for every compact subset of M 
that shares with M a part of its boundary. Namely, let F be a subset of M 
such that V is compact and let F = dV H dM be the part of the boundary 
that V shares with M, then our main result can be stated as follows: 

Theorem. Let (M,g), Rg, hg, V , T, K, c be described as above. Let u > 
be a classical solution of ( tj.g|) . Suppose T is umbilic and there exists eo > 
such that for every point p €zT, g is conformal to the Euclidean metric 6 in 
B(p,eo) r\ M . Assume in addition that for some A > 0, 



(1.2) 




- < K{x) < A, M 
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max ^||VE:(x)||i^(M), ||c||ci(9M)^ < A. 

Then for some C{AI, g, A, y, eg) > and integer m{M, g, A, VjCq) > 0, there 
exist local maximum points of u, denoted as S := {Pi, ...,Pfn}, such that 

distg{Pi,Pj) > ^, ^u{p) < u{Pj) < Cu{Pi), yi / j, 

^ m m 

1=1 1=1 



where 



As a consequence, we have the following energy estimate: 
Corollary 1.1. Under the assumptions of Theorem^ we have 

(1.4) / {\Vgu\^ + u^) < C{M, g, A, V, eo). 

Jv 

If {M,g) is compact. Theorem „ concludes that if the metric is locally 
conformally flat near dM, which is umbilic, then there are only finite blowup 
point in M and the solutions of (jl.2p can be estimated by (jl.Sp . 

The results of Theoremoand the Corollary [TTT] are closely related to a priori 
estimates of the solutions. When {M,g) is compact with a boundary and K 
and c are constants, various compactness results and related discussions can 
be found in ^20j, ^22j, ^30j and the reference therein. The assump- 
tions on K and c in Theorem„should be sharp. The reason is even for the 
interior equations the assumption on K is necessary to have a descrip- 
tion of the blowup phenomenon by (jl.Sp . So Theorem„implies that for this 
three dimensional situation, the assumption on c does not lead to more 
restrictions on K. For dimensions 4 and higher, we expect the situation to 
be much more subtle, first the assumption on K will be more delicate and 
should be made in neighborhoods of its critical points, then the sign of c at 
the blowup point will make a more significant difference. Also, the flatness 
assumptions in the neighborhood of the critical points of K and c might 
be related. As the dimension grows higher, the relations will become more 
subtle, new phenomenon will come out. 

The major step in the proof of„is to establish a Harnack type inequality 
near T. This is done by the well known method of moving planes. Unlike 
most of the previous works on this type of Harnack inequality established 
only for the interior equations, the equation in Theorem„has new features 
that were not dealt with before. Generally speaking, the method of mov- 
ing planes requires a delicate construction of some test functions. Now this 
becomes more difficult because of two reasons, first, the symmetry of the 
domain is destroyed, second, when blowup happens on the boundary or close 
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to the boundary, the function c creates more error terms for the test func- 
tions to control. In this article we found a way to handle all the difficulties 
for dimension three, our approach is motivated by Caffarelli's magnificent 
ideas in [9j [10] on free boundary problems. 

In the next few sections, we mainly focus on the proof of the Harnack 
type inequality on upper half Euclidean balls and we shall mention how to 
obtain the interior estimate away from the boundary T. 

2. Harnack type inequalities on the boundary 

The proof of TheorerOnCan be divided into two parts. First we consider the 
region close to T: N{T, eq). i.e. the eo neighborhood of T. In this region we 
shall establish a Harnack type inequality. Then outside this region, we use 
the techniques developed in ^25j to get the same type of Harnack inequality. 
Note that for this interior part we don't need to assume the metric to be 
locally conformally flat. 

Since (M, g) is locally conformally flat near T, at each point p £ T we 
can flnd a positive smooth function (p such that 6 = (j)'^g where 5 is the 
Euclidean metric. The umbilicity of F implies that dM near p is either a 
piece of sphere or a hyperplane in the new metric. Since both neighborhoods 
are conformally equivalent through the inversion map f{x) = x\x\~'^ , we just 
assume that near p the boundary is a hyperplane, then the equation in the 
neighborhood of p can be written as 

^(f) + ^(f)' = B+CIW' 

dUl) = -c{l? on 

where ei > is a positive constant that only depends on eo, M, 5, F. 
is the upper half ball of radius ei, 9' 5^ = dBf^ n 5M!]_. The computation 
above is based on the following conformal covariant properties of the two 
operators: 

A(-) = 4)-^{^gU-c{n)Rgu) 

u 1 

-9xA^) = ct)~^{dugU + -hgu). 

So the whole equation is reduced to the Euclidean case. From now on in 
this section we just consider the following case (by abusing the notations we 
still use K to denote a positive function on a upper half disk B^^ and c 
to denote a function defined d'B^^^, the lower part of dB^j^): 

on d'B+. 



(2.1) 



f Au + K{x)u^ = 

- — = c(x )u 
^ 0x3 
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where d'B^^ = dB^^ n 5M']_, x' is the projection of x to & B^^. The main 
result of this section is 



Proposition 2.1. Let u G C"^(-B^) n C'^{B^j^) be a positive function that 
solves Ii2.1\) . Assume for some A > 0, K and c satisfy 

(2.2) -^<i^(x)<A, \VK{x)\<AR-\yxe^^ 
and 

(2.3) |c(x')|<A, \Vc{x')\ < AR-\ Vx € S'S+j 
respectively, then for some C(A) > 

(2.4) niaxu • min-u, < C(A)i?^^ 

This Harnack type inequahty reveals important information about the 
interaction of bubbles (large local maximum points) of u which can be seen 
in the following 

Corollary 2.1. Under the assumptions in Theorem \2.1\ there exist a posi- 
tive constant C{A), a positive integer m{A), and a set of m local maximum 
points of u in B^, denoted as S := {Pi, ..,Pm}, such that 

dist{Pi, Pj) > R/C, ^u{Pj) < u{P,) < CuiPj), Vi / j, 

^ m m 



where 



C 

1=1 1=1 



Consequently u satisfies 

(2.6) / (|Vn|2 + n6) < C(A). 



In the proof of Proposition we shall omit a selection process of Schoen 
since it is well known to experts. We only include it in the appendix. At the 
end of this section we indicate the outline of the proof of Corollary 1 2 . II based 
on Proposition 12.11 The proof can be found in ^24j with obvious changes, 
but we indicate the key points for the convenience of the readers. 

2.1. The Proof of Proposition 12.11 We only need to consider the case 
R = \. The general case can be reduced to the case i? = 1 by considering 
the function R^u{R-). The proof is by a contradiction. Suppose ()2.4p does 
not hold, there exists a sequence Ui{xi) such that 

(2.7) Mj(xj) rniniij > i 

B+ 
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where Ui{xi) = max^+ Ui and Xi & Bi. Clearly the above inequality implies 

Ui{xi) oo. By a standard selection process of Schoen [31] and the clas- 
sification theorems of Caffarelli-Gidas-Spruck [8] and Li-Zhu [27J, one can 
consider local maximum of Ui, moreover the following sequence of 

functions 

(2.8) Vi{y) := Ui{xi)~^Ui{ui{xi)~'^y + Xi) 

converge in over any finite domain in the following two cases: 

(1) If \\TCLi^oou'i{xi)xi'i oo, Vi{y) converge in Cf^^i^^) to U{y) which 
satisfies 

/ ^^(y) + 1™^— Ki{xi)U{yf = 0, M^, 
^ ' \ ?7(0) = 1 = maxR3 U. 

(2) If \\Tai^aouf{x.i)Xi^ is bounded, by passing to a subsequence we as- 
sume 

lim uf{xi)xi3 = to. 

In this case Vi converge uniformly to U which satisfies 

r AU{y) + limi^oo Ki{xi)U{yf = 0, ys > ^o, 
(2.10) <^ 93C^(y) =limi^ooQ(x'i)f/3, ys = -^o, 

[ C/(0) = 1 = max U. 

Since the selection process and the application of the classification theorems 
are standard. We put the details in the Appendix. Similar techniques can 
also be found in [TT], [23], [37], etc. 

The proof of Proposition 12. II that follows can be divided into three steps. 
First we rule out case one. i.e. We shall show that the blowup points can 
not be far away from the boundary. In the second step we prove the case of 
limj_+oo c{x',i) < 0. Then in the final step we prove the case of limj^oo c{x^) > 
0. 



2.1.1. Step one. In this subsection we derive a contradiction to 

lim Ti := lim uf{xi)xi3 oo. 

With no loss of generality we assume limj^oo Ki{xi) = 3, in this case the U 
in (|2.9p is of the form 

C/(y) = (l + |y|2)-i 

Let Mi = Ui{xi). In this subsection we use Vi in stead of Vi. i.e. 

Vi{y)=M-^u,{M-^y + Xi). 

The selection process implies B{xi,^) n MS]_ C i?^ so the rescaled domain 
for Vi will have a part of the boundary (the upper part) whose distance to 
is comparable to Mf . By (|2.7p we know on this part of the boundary 
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'^i{y)\y\ > By choosing ^ slowly we can still obtain Vi{y)\y\ > \/i for 
\y\= eiMf. So we define 

and 

= d^i n {yg = -Ti}, d"ni = dUi \ d'Ui. 

Then Vi satisfies 

' Avi{y) + Ki{Mr^y + Xi)vi{y)^ = y e 

(2.11) 

, dsviiy) < Avf yed^in {2/3 = -Ti}, Ti := Mfxis 

Note that in this case we don't need the specific equation for Vi on &Ui, 
By the discussion on the behavior of Vi on we have 

(2.12) Viiy)\y\^oo for y G 
Let A G [5,2] and 

then direct computation gives 

(2.13) Av,\y) + K,{MrY + Xi)vHyf = for |y| > A/2. 

Here we use y'** = A^y/|yp. 

In this step we assume A G [^,2]. Let Hx = Qi\ Bx- Note that for 
simplicity we shall omit i in some notations. Let wx = Vi — vf and we 
consider the equation for wy. 

Awx + bxwx = Qx, Sa, 

(2.14) < 

^ d3wx < A^iwx + A{v^f - dsiv^), &Q.i \ Bx. 
where 6a and are obtained from mean value theorem: 



hx{y) 



5Ki{Mr^y + Xi) [ {tvi{y) + (1 - t)vl{y)fdt 
Jo 

= 3 / {tvi 
Jo 



^i{y) = 3 / {tVi + (1 - t)vtydt. 
and 

QA(y) = (K^iMrY + x.) - K,{Mr^y + x^)){v^{y)f 
is estimated as follows: 

(2.15) \Qx{y)\<Co{A)M-^\y\-\ Sa- 

Since Vi converges in norm over U over any fixed finite domain, is 

close to U^, the Kelvin transformation of U. By direct computation U > 
for \y\ > A and A G (0, 1). On the other hand U < U'^ for A > 1 and \y\ > A. 
So the strategy of the proof is to find a test function hx (i omitted in this 
notation) so that the moving sphere method works for wx + hx, and the hx 
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is just a perturbation of wx, which means h\{y) = o(l)|y|^-'^ in S^. Then 
it is possible to move the spheres from a position less than 1 to a position 
larger than 1 keeping wx + hx > in T,x- But this is a contradiction since 
Wx + hx converges to U — in finite domains and U < U'^ for A > 1 and 

\y\ > A. 

The test function in this section is 

(2.16) hxiy) = CoMr2(_i|y|-2 + A-i|y|-i - Ia-^). 

hx is a radial function, the function hx{r) satisfies 

/ Kir) + Ih'^ir) = -CoMr\-\ r > A, 
\ hx{X) = h',{X) = 0. 

By maximum principle, hx{r) < for r > A. Note that the Cq in (|2.16p is 
the same as the one in (|2.15p . 

Now we consider the equation for wx + hx, from ()2.14|) we have 

(2.17) {A + bx){wx + hx) = Qx + ^hx + bxhx<0 S^. 

The reason is Ahx + Qx ^ and hx < in T,x- 
On d'^li \ Bx we have 

(2.18) d^{wx + hx) < AUwx + hx) + Mv^)^ - ds{v^) + d^hx - A^ihx. 
By the definition of hx one can verify that hx = o[l)\y\~^ over S^. Let 

Ox = {y€ Sa; v^{y) < 2v^{y) }. 

Then we observe that the maximum principle in the moving sphere process 
only needs to be applied over Ox because outside Ox-, wx + hx> 0. Now we 
claim that 

(2.19) K{vff -d^{vl) + d:,hx-K^ihx<^, on ^UnO^. 
Once we have (I2.19|) . (12.180 becomes 

(2.20) d^{wx + hx)<Kii{wx + hx) on &^^r^O~x, 

which is the form for the application of the moving sphere method. 
To see (I2.19p . we first observe that 

3 

dM) = \\y\-^TiVi{y^) + 2\'\y\-''Y,dMy^)yjTi 

i=i 

+\M~^dMy'') on &ni\Bx. 

The second term and the third term on the right are negligible comparing to 
the first term on the right. The reason is Vi ^ U in Cf^^i^), which implies 

Vv,{y^) ^ V[/(0) = 0. 

On the other hand, by assumption Tj — > cxo, we see the first term dominates 
the other two terms. 
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Next we see that the assumption Ti —>■ oo makes |93(Vj^)| dominate A{v^)^ 
as the later is of the order 0{\y\~^). By the definition of hx and the definition 
of the domain Qi, we see easily that 153(^^)1 also dominates d^hx. Finally in 
Ox, is comparable to \y\~'^, so ^ihx is comparable to 0{Mf'^\y\~'^), which 
is much smaller than \d-i{v^)\, (j2.19p is proved. 

The process of making the moving sphere process start is standard, even 
though the boundary condition makes it different from the interior case. 

Lemma 2.1. For any fixed Xq G (^,1) and all large i, 

wxo + /iAo > in Saq. 

Proof of Lemma \2A\ Since Vi —>■ U in Cf^^{R^) and U > for \y\ > A if 
A G (^j !)• By the convergence of Vi to U and the fact hx = o(l) over finite 
domains we can check easily that for any fixed Ri >> 1, 

wxo + ^Ao > in n Br^ . 

Since Ri is chosen sufficiently large, one can find o"(Ao) > such that 

vf°{y)<{l-3a)\y\-^ \y\ > Ri 

and 

> (1 - cr)|yr^ on \y\=Ri. 

The definitions of hx^ and fij imply |/iAo(y)l ^ f 1^1^^ Ivl > -^i- S° *° 
finish the proof of this Lemma it is enough to show 

(2.21) Viiy) > (l - 2a)\y\-^ in O^^- 

On \y\ = Ri and d"Qi we certainly have Vi > {1 — 2a)\y\~^. To apply the 
maximum principle to the super harmonic function Vi — {1 — 2a)\y\~^, we 
need 

(2.22) d3{v,-{l-2a)\y\-^)<A{vf-{l-2af\y\~^). d'Qi 
Using Tj ^ oo we obtain by direct computation that 

93((1 - 2a)\y\-') > A((l - 2a)\yr'f d'n,. 

(j2.22p follows immediately. (j2.22p means the non-positive minimum can not 
be attained on d'Qi. Lemma |2. H is established. □ 



Lemma 12.11 means the moving sphere process can start at A = Aq • The 
purpose of the moving spheres is to show that 

(2.23) wx + hx>Om Sa for all A G [Aq, Ai]. 

Once this is proved, we have U > U'^^ for some \y\ > Ai, which is a con- 
tradiction to the fact that Ai > 1. To see why (|2.23p holds, first, (j2.17p 
and (j2.20p means the maximum principle holds for wx + hx, second, (j2.12p 
means there is no touch on d"Qi. Step one is established. □ 
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2.1.2. Step two. In this step, we deal with the case hnij^oo Ci{x'j) < 0. Let 

—Co = hm Ci{x'j) < 0. 

In this case the hmit of Vi is U, which satisfies 
I' A TT I QrrS 



AC/ + 3C/5 = 1/ G ys > - hm, 



T- 



hni^ — Tj^ , 



U{0) = 1 = max U. 
By Li-Zhu's classification theorem [27], limj^oo Ti = 0. So we define Vi as 

Vi{y) = M-\i{M-^y + x'i) 

where x'^ is the projection of Xi on dM.^. Then f j — > ?7 in Cf^^{M.^). Since 
?7(0) = 1, we know, from Li-Zhu's theorem, that 



(2.24) U{y) = (1 + c^)2 (1 + (1 + 4y\y'\' + (1 + c^)(y3 + 
Let v^{y) = {■^)vi{y'^). In this step we let 



Co ^2^ 



1 + 4' 



Ao = 2(l + Co) 



i Ai = 2(l + cg)-^ 



and we require A G [Ao,Ai] and the moving sphere method will be applied 
on Sa = B~^{0,eiM^) \ B\. Direct computation gives 

r ^vl{y) + K,{Mr^y^ + x'^){vl{y)f = S^, 
(2.25) I 

I dzvf{y) = CiiMr^y^ + x[){vf{y))\ d^x D dRl 

Also by direct computation one can verify that 

' U{y) > U\y), |y| > A, < A < (1 + cg)-^, 

. U{y)<U\y), \y\>\ \>{l + clrl 

For this reason we require A G [Ao,Ai] in this step. The equation for w\ = 
Vi — v} becomes 



(2.26) 



agU-A = Ci(Af-2 . +x'i)i2Wx + Q2, d^x H 

where and ,^2 are obtained by mean value theorem: 



6 



(tvi + (1 - t)v^)Ut, e 



{tvi + (1 - t)vf fdt 



and Qi and Q2 are error terms to be controlled by test functions: 
Qi(y) = {Ki{Mr^y'' + x[)-Ki{M-^y + x',))vf{y)\ Sa, 
Q2(y) = {ci{Mr^y + x'i)-c^{Mr^y^ + x'i))v^{y)\ ^Sa n 
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For Qi and Q2 we use 

\Qi{y)\ < C2WMr^\y\'\ Sa, 

(2.27) |g2(y)| < C2{A)Mr^\y\-\l--^^), ^EaHSM^. 

Note that in Q2 it is important to have the 1 — term, even though we 
don't need this term to appear in the estimate of Qi. 

The construction of the test function here consists of two parts, hi and 
/i2- We first define 

(2.28) hi{r) = C3(A)(r-i - ^A^r'^ - ^A-^)Mr2 
where Cs{A) >> C2 is to be determined. hi{r) is the solution of 

r h'lir) + ^,h[ir) = ^%M-Vi, r > A, 
\ /11(A) = /i; (A) = 0. 

By maximum principle, /ii(r) < for r > A. /ii(|y|) is the first part of 
the test function h\. To define the second part of the test function, we let 
(f) : [1,00) — [0,00) be a smooth non-negative function that satisfies 

r 0(1) = 0, 0(r) = for r > 2, 
\ 0'(1) > 0, 0(r) > for r > 1. 

Let (t)\{-) = 0(A-) and we define 

(2.29) h2{y) = -C4{A)M-^y3M\y\) 
where C4(A) is chosen so large that 

(2.30) C^Mlvl) > 2C2|2/|-=^(1 - ^), \y\ > A. 

\y\ 

Note that on S'Sa := 9Sa n dM.1, 

(2.31) d3h2(.y) = -C^Mr^Mlvl), dshi{y) = d'^x- 
Let hx = hi + h2, the equation for iw;^ + hx is 

(A + Ki{M-^ ■ +x',)Ci){wx + hx) = Qi + Ahx + KiiMf^ ■ +x\)^ihx ^x, 

{83 - Ci{Mr^ ■ +x'i)^2){wx + hx) = -Ci{Mr^ ■ +x'i)^2hx + ^3/12 + Q2 d'^x- 

For the application of the moving sphere method we show that the right 
hand sides of the above are non-positive. Namely we shall show 

(2.32) Qi + Ahi + Ah2 + Ki{Mr^ ■ +x'i)^ihx < 
and 

(2.33) - Ci(Mr2 . +x'i)^2hx + 93/?.2 + <92 < d'J^x n 

where Ox = {y G S;^; Vi{y) < 2v^{y) } is the only place where the maximum 
principle needs to hold, because by the definition of hx, hx{y) = o[l)\y\~^ 
in T,Xj which means wx + hx > in S;^ \ O;^. 
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To see ([232]), first by ([2:2711 and (12:281) one sees that if C3 is large enough 

A/li + Qi < -Jc3M-2|y|-2-5, Sa. 
o 

For /i2 we have 

A/12 = -C4M-2(c/>';((|y|) + ±^'^(\y\))y^ S^. 

So by choosing C3 larger if necessary we have 

A(/li + /l2) + Ql < -^C3M-2|y|-2-5 S^. 

The term KiiM''^ ■ +Xi)^ihx is non-positive. So ()2.32p is proved. 
To see (I2.33p . first by the definition of /12 we have 

93/12 + Q2< -^M-^M\y\), d'J^x. 

Since limj_^oo Ci{x[) = — cq < 0, one uses the fact \y\ < eiMf and the uniform 
bound of the Ci norm of Cj to get 

< o(l)e2|/ll + /l2| <o(l)Mr2|y|-2(l_^) a'SAflOA. 

The second inequality above is because in 0\, ^2 ~ We also use 

h\{X) = in the above. Then (j2.33p follows immediately. 

Next we show that the moving sphere process can get started, namely we 
have 

Lemma 2.2. + /iaq > m Sa,, for all large i. 

Proof of Lemma 12. 2t The proof is similar to the one in step one. For 
\y\ < Ri {Ri » 1)) + /lAo > is guaranteed by the expressions of U, 
[7^0 and the convergence of to U over finite domains. Also hx^ = o(l)|y| 
means Hxq is just a perturbation over finite domains. Moreover on \y\ = Ri, 
there is eo > such that 

U{y)-U^'^{y)>5eo\y-e3\-\ \y\ = Ri, 63 = (0,0,1). 

Also we have 

v^'iy) < (1 - 4eo)(l + cg)"^!?/ - e3\-\ \y\ > Ri,y e Saq 

and 

\U^"{y)-vt'{y)\<eo\y-e3\-\ \y\ > Ri 
for all i large. So for y G Sa \ Br^ we only need to show 

Vi{y) > (1 - eo)(l + Co)~^|y - 631"^ y G Sa \ fi/ji- 
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Since Vi{y) — {l—eo){l+CQ)~^\y—es\~^ is super-harmonic function in Sa\-B_Ri 
and is positive on \y\ = Ri and \y\ = eiMf. The only place to consider is 
d'T,XQ \ Br-^, on which we have 



d3{vi{y)-il-eo)il + cl)~^\y-e3\-^) 
c,{M-^y + x[)vf{y)-i 
Since limj_>oo Ci{x'j) < 0, (j2.34p means 



(2.34) = c,(M-2y + x>f(y)-(l-eo)(l + cg)-i|y-e3rl 



d3ivi{y) - (1 - eo)(l + c^y^y - 63!"') < 

at places where Vi is close to (1 — eo)(l + CQ)~^\y — e^l"^ on \ 
This means 

Vi> {l-€o){l + cl)-^\y-e3\-^ on O'I^x.XBr^. 

Lemma 12.21 is established. □ 

Once we have Lemma [2.21 the moving sphere process can start at A = Aq. 
Note that the equation for wx + hx becomes 

{A + K,{Mr^-+x'^O{wx + hx)<0, S^, 

(93 - Q(Mr2 . +x'^0{wx + hx) < 0, &^x n Ox- 

This means the maximum principle always holds for wx + hx as long as 
it is positive on the boundary d'Tli, which is certainly the case because 
Vi » vf on d'Q.i and hx is a perturbation. So the spheres can be moved to 
Ai = 2(1 + Co)~ 2 to get a contradiction. 

2.1.3. Step three. In this subsection we deal with the case limj^oo ^(2;^) = 
Co > 0. Vi converges uniformly to U in all finite subsets of {y G M^; yz > 
— liTHi^oo Mfxi3 }• U satisfies (still assuming \im.i^ao Ki{xi) = 3) 

/\U + 3^7^ = 0, y3 > - lim.^oo Mfxi^ 
dsU = coU^, y3 = - linii^oo Mfxi^, 
U{0) = 1. 

By Li-Zhu's classification theorem limj^oo M^^Xi^ = cq, and 

[/(y) = (l + |y|2)-i 

Let 

Vi{y) = Mr\i{Mr^y + x'i). 

Then 

Vi Ui{y) = {l + \y - Coe3p)"2. 
Let v^ be the Kelvin transformation of Vi. In this case we let 

Ao = ^(l + cg)^ Ai = 2(l + c2)l 
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because 



Ui>U^ for AG (0,(l + c2)5), |y| > A, 
Ui < for A > (1 + cg)5), \y\ > A. 



Let w\ = Vi — v^, then the equation for w\ is still described by (j2.26p with 
^i! ^2, Q11Q2 defined as in step two and the estimates for Qi and Q2 are still 
(j2.27p . Since \\m.i^oo Ci{x[) = cq > 0, the construction of the test function 
hx is this case is much more delicate. 

The construction of hx consists of two parts. First we shall construct hi to 
control the region close to dBx, then this function is connected smoothly to 
0. Clearly hi creates new difficulties at the regions where it is connected to 
0. Next we use /12 to control the region far enough to dBx- /12 is in regions 
close to dBx and becomes negative as y is far away from dBx- A parameter 
of /i2 is chosen to be large so that /12 not only controls the difficulties from 
Qi and Q25 but also those from hi. One delicate thing is that /12 does not 
create new difficulties. By choosing all the parameters carefully we shall 
obtain the following properties for hy- 



Ahx + KiiM-^y + x'^^ihx + Qi<0 Sa 

dshx - ciiMr^y + x^^/iA + Q2 < 0, d^x n dwl n Ol. 

hx = on dBxHRl, 
(2.35) hx = o{l)\y\-' Sa 



where A > 5Ai is to be determined. Let J7i C M.^ be the region where (j) is 
positive. Let Ti be the level surface of (f). Direct computation shows that 



where Cg is a large constant to be determined. From here we see /12 is C 
across Pi. By direct computation 



We first construct /i2. Let 



4>{y) = 1 - Ai\y + res\ 4, r = \y 



Rl \ B{0, 2A) C ill C \ B{0, .4/3). 



Let 




Rl \ {Qi U Bx) 

T.xnni 



(2.36) 



Ah2 



< 
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From here we observe that /12 is weakly super-harmonic in Y^x since V/12 = 
on Fi. On d^l^i n dM.^ n Ox we need to estimate 

Since 



M 
iTr 

for some M > depending on cq, it is enough to show 

M 
1 + r 

By the definition of /12, if A{cq) is large, we have 



Q(Mr2y + x'i)^2 < -TTZJi Ox n 



(2.37) 93/12 - ^^/i2 + Q2 < 0, n dR% n Oa- 



53/12 - T^/i2 = -2C6Mr2(93 ' ^ 



l + r2 ^ " * ' •''^ 2(1 + r2) 



< -2C6M-2(a3<A-^^Y^)0 
= -2CfiMr2(i^i(2r2)-|r ^ 



2(1 + r2)- 

(2.38) < -—M-'^Air--i(t) OJli n ^ 



Here ^4 is fixed. We observe that if \y\ > 16A, (f>{y) > i. So by choosing Ce 

large enough we have 

(2.39) 

{d3-Ci{M-^-+x'i)C2)h2 + Q2 < -^M-2A3r-t 9(SA\Bi6A)n9M^nOl. 
By (j2.38p we also know 

(2.40) (^3 - Q(Mr2 . +x'i)^2)h2 < d{Yx n B^qa) n SM^. n O^. 

So the most important feature of /i2 is it does not create new difficulties, 
even though it does not control all the error terms in the whole S^. 

Now we define hi to be 
(2.41) 

' -C5M-\r - X)e^y^ YxDBieA 
hi{y) = I Sa\Si7a, 
smooth connection such that /ii < in Sa Pi BnA \ 

where C5 and are large positive constants to be determined. A is the 
one in the definition of /i2, which has been determined. Also the "smooth 
connection" in the definition means |V/ii(y)|, |V2/ii(y)| < Ci{C^, A)M~'^ 
for IQA < \y\ < 17A, y e Sa- 
By the definition of hi, 

Ahiiy) = -C,M-\''y^- + —y3 + N\r-X)) 

r r 

(2.42) < -2C75Mr2r-i yG^isAnSA. 
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Here we require C5 to satisfy 

(2.43) A/ii + Qi < y G n S^. 
On the boundary we have 

d^hi - c^iMr^ + x[)C2hi 
= -C5Mr\r-X){N -CiiUr^y' + x'i)C2) ^(Sa n Sisa) n ^R^ n O^. 
So by choosing N{cq) large and C5 larger if necessary we have 

(2.44) d^hi - a {Mr^y' + x'J^a/ii + Q2 < d{^x n Sisa) n n O^. 

So far CsjA^, A are all determined, we finally determine Cq. One last 
requirement for Cq is to control the bad part of hi in {Bua \ Biqa) H Sa- By 
choosing Ce larger if necessary we see from (12.38P and (12.39P that the errors 
caused by bending hi to can be controlled by /12. Therefore the first two 
equations of (|2.35p are established. Here we also use the fact that hx < 
in Tix. By the definitions of hi and h2, the other two equations of (j2.35p 
are also satisfied. We are left with how to let the moving sphere process 
start. This part is similar to step two. We leave the details to the interested 
readers. Proposition 12 . II is established. □ 

Proposition 12 . II is the major step in deriving the energy estimate in Corol- 
lary 12. 1[ Based on the previous works of Y. Y. Li p3] , Han-Li [22] and Li- 
Zhang [24j . it is standard to derive Corollary [27TJ Therefore we only mention 
the major steps and the main idea in this argument: First one uses Schoen's 
selection process to find all large maximums of u in, say, -B^. Around each 
of these local maximums there is a small neighborhood in which u looks 
like a standard bubble which has most of its energy in it. The distance 
between these local maximums is the crucial information to find and this 
is the place where the Harnack inequality in Proposition 12.11 is used. The 
essential difference between this locally defined equation and those globally 
defined equations (such as the ones in [23]) is that one can not find two 
bubbles closest to each other. For each bubble, there certainly exists a bub- 
ble closest to it, but one certainly can not assume the first bubble is the 
closest one to the second bubble. This difficulty, which comes from the lo- 
cal nature of the equation, requires a different approach than those in [233 • 
The way to overcome this difficulty is to rescale the equation so that after 
the scaling, the equation is centered at the first bubble and the distance 
between the first bubble and the second bubble is one. Then the Harnack 
inequality in Proposition 12.11 applied to the local region implies that the first 
bubble and the second bubble must have comparable magnitudes. Then it 
is possible to show that two bubbles can not tend to the same blowup point 
because otherwise a harmonic function with a positive second order term 
can be found. This second term will lead to a contradiction in the Pohozaev 
Identity. Note that the important second order term can only be proved 
to be positive if the adjacent bubbles have comparable magnitudes, which 
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is the key information revealed by the Harnack inequaUty. Once we have 
known that ah bubbles are far apart, it is possible to use standard elliptic 
estimates to show that the behavior of u near each large local maximum is 
like a harmonic function with fast decay. So (12.51) as well as (|2.6p can be 
obtained. □ 

2.2. Harnack inequality in the interior of M. Now we consider V \ 
A^(r,eo). Let d be the distance between V \ NiV^eo) and dM. Over this 
interior region we don't need to assume M to be locally conformally flat. 
We shall establish the following inequality: 

Proposition 2.2. Let p G V\N{T,eQ), there exist e2{M,g,A,d) > and 
C{M, g, A,d) > such that 

(2.45) sup u- inf n < ^ 

B(p,t2/3) -B(p,2e2/3) 

The proof of Proposition 12.21 follows from the argument in [25j . Only 
small modification is needed to adjust to the current situation. The outline 
is as follows: Suppose (j2.45p does not hold, then we can select Xi as a 
sequence of local maximums of Uj. Then the equation can be written in a 
conformal normal coordinates centered at Xi. After rescaling, Ui becomes Vi, 
a sequence that converges in Cf^^{M.^) to a standard bubble whose maximum 
is 1. Then consider the Kelvin transformation of fi(y): v^{y). By comparing 
the equation for Vi and we see that the only new term in this context is 
the term: 

{K^{M-^y^ + Xi)-Ki{M-^y + x^)){v^f in 

where Ea becomes a symmetric domain defined appropriately. The above 
term is of the harmless order 0(M^~^|y|~^). By using the same test function 
we used in [25], a contradiction can be obtained correspondingly. We leave 
the details to the interested readers. Proposition 12.21 is proved. □ 

Proposition 12.11 and Proposition 12.21 lead to Theorem„by the arguments 
in [24] and [26]. 

3. Appendix 

In this section we provide some details of the Schoen's selection process 
for the convenience of readers. We only provide details for case one and case 
two. The details for case three are similar. 

3.1. A Calculus Lemma. First we use a Calculus lemma to simplify some 
computations. This Calculus lemma was first used in [24]. 

Lemma 3.1. Let u G C°{Bi n {t > -T}) T > 0. Then for every a > 0, 
there exists x £ Bi n {t > —T} such that 

1 1 

u(x) > — max u, and a'^u(x) > — u(0). 



18 LEI ZHANG 

where = B{x,a) n {t > -T], a = {I - \x\)/2. 

Proof of Lemma I3.lt First we remark that if T > 1, the selection is over 
the whole Bi. To prove the lemma, consider 

v{y) = {l-\y\Tu{y). 

Let X G Bi be a maximum point of v and let o" = (1 — |x|)/2. By comparing 
v{x) with t>(0) and v{z) for all z in BJ{x), we see that x and a have the 
desired properties. □ 

3.2. The selection process for case one in the proof of Proposition 

I2.lt In this subsection we explain why be considered local max- 

imum of Ui, assuming n?(xj)xj3 oo. We shall apply Lemma [3TT] for T > 1. 

1 

i.e. the selection is over the whole Bi. Let = |xj3|/2, so Ui{xi)rf — > oo. 
Apply Lemma |3. II to the function Ui{xi + ri-) over Bi, then the conclusion 
can be translated as follows: there is Oj G B{xi,ri) such that 



1 



Ui{ai) > max 2 2Uj(x) 

B{ai,ai) 

where ai := ^(r^ — \ai — Xi\) and 

(ylui{ai) > 2 ■iUi{xi)rf oo 

Let 

Vi{y) = Ui{aiy^Ui{ui{aiy^y + a^) < Ui{aifai 

Then 

Vi{y)<22 for |y| < Ui(aj)^cri ^ cx) 

Clearly satisfies 

A^idz) + Ki{ui{ai)~'^y + ai)vi{yf = |y| < Ui{aifai 

So by standard elliptic theory we know there is a subsequence of Ui(y) (still 
denoted by Vi{y)) that converges uniformly to Uo{y) on all compact subsets 
of M'^. Uo{y) satisfies 

(3.1) AUo{y) + KoUoiyf = V^R^ 

where Kq = limj__^oo Ki[ai). With no loss of generality we assume Kq = 3, 
so by the well known classification theorem of Caffarelli-Gidas-Spruck [8], 

= ( i ^ 2^ 

1 + — 

for some z G and /U > 1. So Uo{y) has an absolute maximum at z. 
Consequently Vi{y) has a local maximum at Zi close to z when i is large. 
Let Xi = Ui{ai)~'^Zi + a^, then {icj} are local maximum points of Ui, also it 
is easy to verify that 

Ui{xi) > ^Ui{ai) > f^u{xi) 
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and Xi G B{ai,ai) C B{xi,ri) which means Xi-^ > ^Xi^, consequently 

uf{Xi)Xi3 oo 



So we can consider at the beginning. The Vi{y) defined in ([2 

-'loc 



converges in C? (M^) to 



C/(y) = (l + |y|^)-2. 



3.3. Selection for case two in the proof of Proposition I2.lt In this 
subsection we show that be considered as a local maximum of it,- 

under the assumption uf{xi)xi3 = 0(1). 

By applying Lemma [3T1 to nj(xj + i-) with a = 1/2 and T = —Xi^i we can 
find Xi £ B{xi, 1/2) n {x-^ > —Xi'^} such that 



Ui{xi) > 2 2 maxuj over B{xi,ai) D {x^ > —Xi^} 



1 



1 



and afui{xi) > 2 ■2Ui(xi) where o"j = 2(2 ~ ~ these two in- 

equalities we have 



1 



Ui{xi) > Ui{xi) 00, Ui{xi)af —>■ 00 

and Ui{xi) infg^+ Ui > i. Now we can assume that % := uf{xi)xi3 is bounded 
because otherwise by replacing Xi by Xi we go back to the case one, which 
has been discussed. Let 

Viiy) = Ui{xi)~^Ui{ui{xi)~'^y + Xi). 

Then Vi satisfies 

Mi{y) + Ki{ui{xi)-'^y + Xi)vi{yf = y G Qi _ 
dzvi{y) = aimixiy^y' + x'i)vi{yf {ya = -Ti} n dQi 

where Cli = {y; Ui{xi)~'^y + G B^}- It is clear that the lower part of 

1 

dQi is a subset of {ys = —Ti}. Moreover by the facts Ui{xi)x^^ < C and 
1 

Ui{xi)a? 00 we know Vi is uniformly bounded on all compact subsets of 
{2/3 > — limi^oo Ti}. With no loss of generality we assume that Ki{xi) — > 3, 
then by elliptic estimates, Vi{y — TiC^) converges in norm on all compact 
subsets of M"^ to C7, which solves 

f Af7 + 3C/^ = 

\ d^tJ = cm yg = 

where c = Imii^^ Ci{x'^) . By Li-Zhu's classification result [27], we know 
that U has a unique maximum point zq in M5]_. By the convergence 
of Vi{- — TiC^) to U, we can find {yi}i=i,2,.. as local maximum points of 
Vi{- — TiCs) that approach zq as i ^ 00. Then by the definition of Vi, we 
know that Ui{xi)~'^{yi — Tie^) + Xi are local maximum points of Ui. So we 
redefine Xi SlS Hi (xi) "^{yi — TjCg) + Xj. So the Vi defined in (j2.8p converges 
to the function Ui in (|2lli]l . 
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